Veze izmedju elemenata
clear variables
syms omega L C w
w = omega
w = 
ω
 
C = 1/(L*w^2)
C = 
Pobuda
syms Um s t
Ug = laplace(Um*sin(w*t), t, s)
Ug = 
Jednacine
syms I0 IL U0 UC
jednacine = [Ug == L*(s*IL - I0) + UC, IL == C*(-U0 + s*UC)]
jednacine = 
 
odziv = solve(jednacine, [UC, IL])
odziv = struct with fields:
UC: (I0*L*omega^4 + Um*omega^3 + I0*L*omega^2*s^2 + U0*omega^2*s + U0*s^3)/(omega^2 + s^2)^2 IL: (I0*L*omega^2*s - U0*omega^2 + Um*omega*s + I0*L*s^3 - U0*s^2)/(L*(omega^2 + s^2)^2)
 
odziv.UC
ans = 
Inverzna Laplasova transformacija
syms uC(t)
 
uC(t) = ilaplace(odziv.UC, s, t)
uC(t) = 
 
vrednosti = [I0 == 0.2, L == 0.1, U0 == 2, Um == 2, w == 2*pi*50]
vrednosti = 
 
uC_num(t) = subs(uC(t), lhs(vrednosti), rhs(vrednosti))
uC_num(t) = 
Crtanje grafika
 
fplot(t, uC_num(t), [0, 0.7])
xlabel('t');
ylabel('uC(t)');
grid on;